NASA Technical Memorandum 100710 


Atmospheric Transmission Coefficients 
for Laser Pulses with Spectral Widths of 
a Few Tenths of a Nanometer, Over the 
Wavelength Region from 800 to 860 
Nanometers 


H.G. Safren 


JULY 1988 


(*ASA-T«-1CG7 10) ATHCSFHEBIC IBANSFISSICN N89-1£1£3 

CCEFFICIES1S f Oi LASEB PULSES II1H SPEC1BAL 
k JD1HS CF A FEK IE M1HS CF A K A AC EET EB CVEB 

Ihl iAVELEBGIE BEG ICS FBCH 8CC 1C 860 Unclas 

AABCHE1BBS (SASA) 18 p CSC1 09 F G3/17 0187240 




NASA Technical Memorandum 100710 


Atmospheric Transmission Coefficients 
for Laser Pulses with Spectral Widths of 
a Few Tenths of a Nanometer, Over the 
Wavelength Region from 800 to 860 
Nanometers 


H.G. Safren 

Goddard Space Flight Center 
Greenbelty Maryland 


NASA 

National Aeronautics and 
Space Administration 
Goddard Space Flight Center 
Greenbelt. Maryland 20771 


1988 



CONTENTS 


Page 


INTRODUCTION 1 

ANALYSIS 2 

DESCRIPTION OF COMPUTER PROGRAM 2 

PLOTS OF EFFECTIVE ATMOSPHERIC TRANSMISSION 3 

SUMMARY AND CONCLUSIONS 9 

REFERENCES 10 

APPENDIX A — COMPUTER PROGRAM LISTINGS 11 


PRECEDING PAGE BLANK NOT FILMED 




111 



INTRODUCTION 


For laser communication experiments between space and ground terminals, the spectral width of the optical 
source can have a significant effect on the atmospheric transmission. The spectral width of a pulsed source 
depends both on the temporal length of the pulse and on the frequency stability of the laser, and also on the 
frequency spread between lasers if more than one laser is used to achieve the desired transmitter power. 

In the 800 to 860-nanometer spectral region, the atmospheric absorption lines (almost entirely water lines) 
have linewidths of a few hundredths of a nanometer. Thus, for a laser pulse with a spectral width of several 
tenths of a nanometer, such as those produced by AlGaAs laser diode arrays, the total atmospheric transmis- 
sion cannot be accurately computed by using the transmission coefficient at a single wavelength. In such 
cases, it becomes necessary to take a weighted average over the pulse spectral width, with the weights de- 
pendent on the pulse spectral shape. The effective transmission coefficient thus computed will depend on the 
pulse’s center wavelength and on its spectral width and shape. 

In this report, the effective atmospheric transmission coefficient is computed as a function of the pulse 
center wavelength over the wavelength region from 800 to 860 nanometers. The computations were per- 
formed for several pulse spectral widths and for a gaussian spectral shape. Computations were also per- 
formed for a spectrally flat pulse to evaluate the sensitivity of the spectra to the pulse shape. It turns out 
that the flat-pulse spectra differ only slightly from the gaussian-pulse spectra, which indicates that the exact 
pulse shape is not critical. 

The single-wavelength (zero spectral width) transmission spectrum over the 800 to 860-nanometer region is com- 
puted in a line-by-line manner. The absorption line catalog used is the 1982 edition of the absorption line 
tape compiled by the Air Force Geophysics Laboratory, which lists all the known absorption lines for the seven 
most important atmospheric absorbing species, from the microwave region to the visible region (References 
1 through 5). The computer program used to perform the calculations is described in Reference 6. 

As would be expected, the effective atmospheric transmission spectrum becomes progressively smoother as 
the spectral width of the pulse increases, so that in general, a spectrally wide signal suffers less energy loss 
as it traverses the atmosphere than does a spectrally narrow signal. Thus, the use of a wider pulse may give 
greater flexibility in choosing the operating frequency of a direct detection laser communication system. 
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ANALYSIS 


Let e(X) denote the energy spectrum of the pulse in the absence of atmospheric attenuation, where X is the 
wavelength. Then, if t(X) denotes the transmission coefficient at the wavelength X, the total received pulse 
energy is given by: 


Received pulse energy = j e(X)t(X) dt. (1) 

PULSE 

If this is equated to the energy that would be received in the absence of atmospheric attenuation multiplied 
by an effective pulse transmission coefficient teff, then the effective coefficient may be written as 


j e(X)t(X) dX 

Wf = (2) 

j e(X) dX 

PULSE 

Clearly, the actual energy e is not important; the energy will therefore be normalized relative to its value at 
the pulse center frequency. Thus, if the quantity s = e(X)/e(AJ, then the above equation may be written as 

j s(X)t(X) dX 

W = ^ (3) 

J s(X) dX 

PULSE 


where s(X) may be regarded as the pulse spectral shape. 

DESCRIPTION OF COMPUTER PROGRAM 

To compute the integrals in Equation 3, it is sufficiently accurate to use a simple trapezoidal rule; thus, the 
effective coefficient may be approximated by: 

» — 2'^(Sntn "F S n + 1 t n + 1 )• (X n + 1 X„) 

Lff s 

Z'M s„ + s n+I )»(X n+I - X n ) 

where s n and t n denote s(X n ) and t(X„), respectively, and the summation is taken over a grid of points X n 
which spans the spectrum of the pulse (and need not be equally spaced). The grid of points is constructed by 
taking an equally spaced grid, with a spacing small enough to cover each absorption line by several points, 
and then adding to the grid the centers of the absorption lines, to ensure that the bottoms of the lines are not 
skipped over. 


Two lineshapes are used in the calculations: a fiat (constant) shape, and a gaussian shape. For the gaussian 
shape, the grid of points extends two halfwidths to either side of the line center. 
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To calculate an effective atmospheric transmission spectrum for a given pulse shape and width, a grid of 
points is constructed as described above, spanning the entire wavelength region from 800 to 860 nanometers, 
by taking an equally spaced grid and intercalating the center wavelengths of all the absorption lines. The 
center wavelength of the pulse is then stepped along this grid; at each point, the effective transmission 
coefficient of the pulse is computed from Equation 4. 

The computed transmission coefficients are stored in a disk file for later plotting. Some representative plots 
are shown in the next section of this report. 

PLOTS OF EFFECTIVE ATMOSPHERIC TRANSMISSION 

Due to the large amount of spectral detail in the 60-nanometer-wide region considered in this report, the 
plots of the effective transmission spectra are divided into three subregions, each 20 nanometers wide. 

Figures 1, 2 and 3 show the effective transmission spectra over these three subregions, for pulses having a 
gaussian spectral shape and spectral widths (full width at half maximum) of 0.1, 0.3 and 0.5 nanometers. In 
each plot, the single-wavelength transmission spectrum (which may be regarded as the spectrum for pulses 
of zero spectral width) is also shown. From these three plots it is seen that the spectrum is quite flattened 
out for pulses which are 0.5 nanometers wide. 

Figure 4 shows a portion of the same spectra, over the smaller region from 815 to 820 nanometers (with an 
expanded scale), to more clearly show the details of the spectra. The computer program is capable of plot- 
ting the spectra over any subregion. 

Figure 5 shows the spectra for 0. 1 , 0.3 and 0.5-nanometer-wide pulses, with both gaussian and flat spectral 
shapes, over the 815 to 820-nanometer region. While the spectra for the two shapes are noticeably different, 
the differences are relatively minor, which indicates that the transmission spectrum is not very sensitive to 
the exact spectral shape of the pulse. 
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PERCENT TRANSMISSION 
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Figure 1. Transmission Spectra for Laser Pulses of Various Spectral Widths— Region 1 
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PERCENT TRANSMISSION 



Figure 2. Transmission Spectra for Laser Pulses of Various Spectral Widths— Region 2 
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PERCENT 




840 845 850 855 860 

WAVELENGTH OF CENTER OF LASER LINE (nm) 

Figure 3. Transmission Spectra for Laser Pulses of Various Spectral Widths — Region 3 


6 





PERCENT TRANSMISSION 
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Figure 4. Transmission Spectra for Laser Pulses of Various Spectral Widths— Expanded Subregion of 
Region 1 
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SUMMARY AND CONCLUSIONS 

Calculations were made to determine the effective transmission spectra, over the wavelength region from 

800 to 860 nanometers, for laser pulses with spectral widths of 0.1, 0.3 and 0.5 nanometers. 

Several conclusions may be drawn from the results of this study: 

1 . For laser pulses with spectral widths of several tenths of a nanometer, the effect of atmospheric ab- 
sorption lines depends on the spectral width of the pulse; 

2. Atmospheric water vapor absorption lines can be very deep, but they are only a few hundredths of a 
nanometer wide in this region of the spectrum, so that laser diode transmitters may not be severely 
affected even if they are centered on a line, or overlap one or more lines; 

3. Pulse spectral widths of several tenths of a nanometer may reduce atmospheric absorption by water 
lines to as little as 10 or 20 percent, even when the laser is centered on a line (the remaining extinc- 
tion of about 50 percent shown in the plots is due mainly to aerosol absorption and scattering); and 

4. The pulse spectral shape appears to have a rather minor effect on the transmission spectrum. (The 
spectra for flat pulses are somewhat more jagged than those for gaussian pulses, which might be ex- 
pected because the flat pulse, as it is moved along the wavelength axis, includes and excludes absorp- 
tion lines much more suddenly than does a gaussian pulse, with its long tails.) 
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APPENDIX A 

COMPUTER PROGRAM LISTING 


PROGRAM SMOOTH 


C 
C 
C 
C 
C 
C 
C 
C 

C***************************************************************************** 

c 

BYTE INFILE(15), OUTFIL(15), VMFILE(15), 

* SHAPE, GAUSS, CONST 


This routine takes a transmission spectrum, in which each value is the 
transmission coefficient at that exact wavelength (in Angstroms, for this 
routine), and replaces that value by the effective transmission coefficient 
for a laser line of finite width (with a constant or gaussian spectral 
shape) centered at the given wavelength. The resulting spectrum is a 
“smoothed” spectrum for the effective transmission coefficient. 


INTEGERS DATCNT 


REAL *4 LOWERW, WAVLEN, TRANS 

DATA INFILE/ 15 *“000/, OUTFIL/15*“000/, VMFILE/15*”000/, 

GAUSS/ 'G 7 , CONST/ 'C 7 


Ask user for the spectrum file to be "smoothed” 


TYPE 100 

100 FORMAT(//t2, Transmission spectrum file to be “smoothed”: ',$) 

ACCEPT 101, INFILE 

101 FORMAT(15Al) 


Ask user for the name of the “smoothed” file 


TYPE 200 

200 FORMAT(//t2, 'Name of file to hold the “smoothed” spectrum: ',$) 

ACCEPT 101, OUTFIL 


Ask user for the lineshape 


TYPE 250 

250 FORM AT(//t2, 'Lineshape — gaussian or constant (G or C)? ',$) 
ACCEPT 251, SHAPE 

251 FORMAT(Al) 
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Ask user for the laser line width 


TYPE 300 

300 FORMAT*// 

* t2, Halfwidth (at halfmax for gaussian shape)'/ 

* t2, 'of laser line, in Angstroms: ',$) 

ACCEPT 301, HWIDTH 

301 FORMAT(F6.2) 


Open the input and output files; 

copy the header information to the output file. 


OPEN (UNIT = 1 , NAME = INFILE, TYPE = 'OLD', DISP='KEEP', 

FORM = 'UNFORMATTED', RECORDSIZE=2) 

OPEN (UNIT =2, NAME =OUTFIL, TYPE = 'NEW ', DISP='KEEP', 

FORM = 'UNFORMATTED', RECORDSIZE = 2, INITIALSIZE = 100) 


READ(l) 

PLTCOD, 

ANGMIC 

READ(l) 

LOWERW, 

UPPERW 

READ(l) 

ATL, 

ATU 

WRITE(2) 

PLTCOD, 

ANGMIC 

WRITE(2) 

LOWERW, 

UPPERW 

WRITE(2) 

ATL, 

ATU 


Open a direct access working file in VM to hold a 
copy of the input file. Copy the rest of the input file 
(the data points) to the VM file; count the number of 
data points. 


OPEN (UNIT =3, NAME='VM: WORK. DAT', TYPE= 'NEW', DISP= 'DELETE', 

* ACCESS = 'DIRECT', FORM = 'UNFORMATTED ', RECORDSIZE=2, 

* INITIALSIZE = 100) 

ICOUNT = 0 

400 READ* LEND =450) WAVLEN, TRANS 
ICOUNT = ICOUNT+ 1 
WRITE(3 ICOUNT) WAVLEN, TRANS 
GO TO 400 

450 DATCNT = ICOUNT 

CLOSE (UNIT = 1 , DISP= 'KEEP') 


Position laser lineshape at high wavelength end of spectrum, 
so that its high-wavelength, two-halfwidth point “hits” the 
upper edge of the data region (recall that the spectrum is 
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stored “backwards” in the files, from high to low wavelengths) 


READ (3'1) WAVLEN, TRANS 
NHIGH = 1 
AHIGH - WAVLEN 
ACENTR WAVLEN - 2. *H WIDTH 

NCENTR - 2 

500 IF (NCENTR .GT. DATCNT) GO TO 5000 
READ (3 'NCENTR) WAVLEN, TRANS 
IF (WAVLEN .LE. ACENTR) GO TO 510 
NCENTR = NCENTR + 1 
GO TO 500 

510 ACENTR = WAVLEN 

ALOW « ACENTR - 2. *H WIDTH 
NLOW » NCENTR + 1 

520 IF (NLOW .GT. DATCNT) GO TO 5000 
READ (3 mOW) WAVLEN, TRANS 
IF (WAVLEN .LE. ALOW) GO TO 530 
NLOW = NLOW + 1 
GO TO 520 

530 ALOW = WAVLEN 


Step laser lineshape along the grid of points; 

at each step, compute the effective transmission coefficient. 


1000 NCENTR = NCENTR + 1 ! Next center point 

READ (3 'NCENTR) WAVLEN, TRANS 
ACENTR « WAVLEN 

1 100 NHIGH - NHIGH + 1 ! Next high point 

READ (3TMHIGH) WAVLEN, TRANS 
IF ((W A VLEN-ACENTR) .LE. 2.*HWIDTH) GO TO 1150 
GO TO 1 100 

1150 NHIGH - NHIGH - 1 

IF (NHIGH .EQ. 1) NHIGH = 2 ! MUST have NHIGH .GE. 2 

READ (3 'NHIGH) WAVLEN, TRANS 
AHIGH = WAVLEN 

1200 NLOW = NLOW + 1 ! Next low point 

IF (NLOW .GT. DATCNT) GO TO 5000 
READ (3 mOW) WAVLEN, TRANS 
IF ((ACENTR-WAVLEN) .LE. 2.*HWIDTH) GO TO 1220 
IF ((ACENTR-WAVLEN) .GT. 2.*HWIDTH) GO TO 1230 

1220 NLOW = NLOW + 1 
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IF (NLOW .GT. DATCNT) GO TO 5000 
READ (3 'NLOW) WAVLEN, TRANS 
IF ((ACENTR- WAVLEN) .LE. 2.*HWIDTH) GO TO 1220 
ALOW = WAVLEN 
GO TO 1300 
C 

1230 NLOW = NLOW - 1 

READ (3 mOW) WAVLEN, TRANS 
IF ((ACENTR-WAVLEN) .GT. 2.*HWIDTH) GO TO 1230 
NLOW = NLOW + 1 
READ (3 mOW) WAVLEN, TRANS 
ALOW = WAVLEN 
C 

1300 SUMNUM = 0. ! Calculate sums 

SUMDEN = 0. 

C 

DO 1350 I = 1 ,NLOW - NHIGH + 1 
II = I 

N = NLOW - (II- 1) 

READ (3 W) WAVLEN, TRANS 

READ (3 'N- 1 ) WVNEXT, TRNEXT 

IF (SHAPE .EQ. GAUSS) GO TO 1310 
IF (SHAPE .EQ. CONST) GO TO 1315 
C 

1310 EXPN = ((WAVLEN-ACENTR)/HWIDTH)**2 * ALOG(0.5) 
EXPN1 = ((WVNEXT-ACENTR)/HWIDTH)**2 * ALOG(0.5) 
SHAPN = EXP (EXPN) 

SHAPN1 = EXP (EXPN I) 

GO TO 1320 
C 

1315 SHAPN =1. 

SHAPN 1 = 1. 

C 

1320 TRMNUM = 0.5 * (SHAPN*TRANS + 

* SHAPN UTRNEXT) * (W VNEXT-W A VLEN) 

TRMDEN = 0.5 * (SHAPN*SHAPNl) * (WVNEXT-W AVLEN) 
SUMNUM = SUMNUM + TRMNUM 

SUMDEN = SUMDEN + TRMDEN 

1350 CONTINUE 
C 

TEFF = SUMNUM/SUMDEN 
WRITE(2) ACENTR, TEFF 
GO TO 1000 


Close the files 


5000 CLOSE (UNIT =2, DISP= 'KEEP') 
CLOSE (UNIT =3, DISP = 'DELETE 0 
STOP 
END 
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